Kinetics of Ring Formation 
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We study reversible polymerization of rings. In this stochastic process, two monomers bond and 
as a consequence, two disjoint rings may merge into a compound ring, or, a single ring may split into 
two fragment rings. This aggregation-fragmentation process exhibits a percolation transition with a 
finite-ring phase in which all rings have microscopic length and a giant-ring phase where macroscopic 
rings account for a finite fraction of the entire mass. Interestingly, while the total mass of the giant 
rings is a deterministic quantity, their total number and their sizes are stochastic quantities. The size 
distribution of the macroscopic rings is universal, although the span of this distribution increases 
with time. Moreover, the average number of giant rings scales logarithmically with system size. 
We introduce a card-shuffling algorithm for efficient simulation of the ring formation process, and 
present numerical verification of the theoretical predictions. 
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I. INTRODUCTION 

Percolation [l|, Q controls many natural processes from 
olymer gelation and diffusion in porous media 

0] , to the spread of forest fires d, Q or infectious dis- 
eases [Toj - Tl^ ] . In the standard percolation picture, a sys- 
tem evolves from a state with many small, microscopic, 
clusters into a state with a single, macroscopic, system- 
spanning, cluster. This phase transition is continuous, 
and it is controlled by the total number of connections 
between elementary units in the system. 

In this study, we show that restricting the structure 
of the clusters leads to a different percolation behavior 
where multiple macroscopic clusters coexist. Percolation 
with multiple giant clusters has been recently reported 
in theoretical studies [HI, [3], and it is relevant to the 
production of colloidal micro-gels [H, HH . 

Our starting point is the classic polymer gelation pro- 
cess introduced by Flory [l?], [HI , a simplified model 
that is essentially the mean-field theory for percolation 
0) Qjrlllj]. In this polymerization process, a very large 
number of molecular units ( "monomers" ) join irreversibly 
to form clusters ( "polymers" ) . This process has a second 
order phase transition between a "sol" phase in which all 
polymers are finite to a "gel" phase in which a single gel 
containing a finite fraction of the monomers in the sys- 
tem emerges. With time, this gel grows and eventually, 
it engulfs the entire system. 

In the Flory model, there is no limit on the num- 
ber of bonds per monomer, and the resulting polymers 
may have very different structures. We modify the poly- 
merization process so that all polymers have the same 
structure. In our version, all monomers have exactly 
two bonds, so that all polymers are rings. Rings oc- 
cur in magnetized powders or beads [24| because due to 
magnetic interactions, linear chains are unstable with re- 
spect to formation of rings (Figure [T]) . As is the case 
for magnetic beads, we consider directed rings where the 
bonds have directionality (Figure The results extend 



FIG. 1: Rings made of magnetic beads. 



to undirected rings. 



II. 



AGGREGATION-FRAGMENTATION 
PROCESS 



At time t = 0, our system consists of N isolated 
monomers. These particles bond to form polymeric rings 
through the following process. In each elementary step, 
two monomers are selected at random and a first bond 
is drawn between them. Subsequently, both monomers 
drop an existing bond and then, the two "dangling" 
monomers form a second bond, as shown in Figure [5] 
Time is updated, t -> t + At with At = 2/N, after 
each step so that each monomer experiences one bond- 
ing event per unit time. We note that the directionality 
of the first bond dictates the directionality of the second 
bond. This polymerization process conserves the total 
number of bonds because two bonds are gained and two 
bonds are lost in each event. We assign an imaginary 
self-bond to every original monomer, so that formally, 
the original monomers have a ring structure. Therefore, 
the total number of bonds in the system equals N. With 
this formulation, the polymerization process maintains a 
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FIG. 2: (a) Inter-ring bonds lead to aggregation. 
ring bonds result in fragmentation. 



(b) Intra- 



ring topology as every monomer has exactly two bonds. 

The above polymerization process is equivalent to an 
aggregation-fragmentation process. When a monomer 
that belongs to ring of size i bonds with a monomer that 
belongs to a different ring of size j, a composite ring 
with size i + j forms (Figure^). Hence, rings undergo 
the aggregation process 



K, 



i+j 



with 



Kij = ij. 



(1) 



The aggregation rate K%j is proportional to the product 
of the sizes because there are i x j distinct pairs that can 
bond. We note that the aggregation process (p} alone 
constitutes the Flory model. 

A bond between two monomers in the same ring breaks 
that ring into two smaller rings. Schematically, the frag- 
mentation process is (Figure [2p) , 



i+j 



1,3 



with 



i + j 
N 



(2) 



Due to the circular symmetry, the fragmentation rate F+j 
is proportional to the ring size, while the factor 1/N re- 
flects that for fragmentation to occur, one must pick two 
monomers within the same ring. Also, the aggregation- 
fragmentation process specified by JT]) and © is re- 
versible because for every aggregation event, there is an 
opposite fragmentation event, and vice verse. 

Let rfc(i) be the density of rings made of k monomers 
at time t. That is, if Rk is the expected number of rings 
of size k, then rk = Rk/N. This size density obeys the 
rate equation 



dt 



J2 

i+j=k 



1 

N 



J r 3 



fc(fe-l) 



(3) 



The first two terms represent gain and loss due to the 
aggregation process (J} , and the last two terms represent 
gain and loss due to the fragmentation process ([2]). In 
particular, let us consider the two loss terms. The to- 
tal aggregation rate is, by definition, the ring size k, but 
the total fragmentation rate, Fk = Y]. , ._ 



k Fij, grows 



quadratically with size, Fk 



±( k ) 



_ fc(fc-l) 



Our 



N\2) ~ 2N 

goal is to understand the time evolution of the density 
rk(t), starting from the monomer-only initial condition, 
rfe(O) = 4,o- 



III. FINITE RINGS 

Our implicit assumption is that the system is very 
large. When N — > oo, one can use perturbation the- 
ory with the inverse system size being the small param- 
eter [2||. We expand the size distribution to first order, 
Tk = Cfe + -kgh + • • • j and substitute this form into ([3]) to 
obtain the rate equation 



dc k _ 1 ^ 
~~dT ~ 2 



IJCiCj 



(4) 



i-\-j—k 



The initial condition is Cfc(O) = Sk,i- The two terms in 
this equation describe gain and loss due to aggregation. 
To zeroth order, the fragmentation process is negligible 
because the likelihood of picking two monomers within 
the same ring vanishes when N — > oo. Equations (JJ) de- 
scribe the evolution of the size distribution in the Flory 
model, where there is no fragmentation. The solution to 
these equations is well-known (see [2(| [26[ for a review) 



Cfc(t) 



1 



k ■ kl 



{ktf 



(5) 



Let M n (t) = Ylk>x k ric k(t) be the nth moment of the 
distribution Ck- The second moment diverges at a fi- 
nite time as M2(t) = (1 — t) -1 for t < 1, a signa- 
ture of the percolation transition at time t = 1. The 
first moment, M\(i), provides additional information 
about this phase transition. Consider the "missing mass" 
g(t) = 1 — Mi(t). This quantity obeys the transcenden- 
tal equation 



1 



-gt 



(6) 



When t < 1, there is only the trivial solution g = 0, 
and hence, finite clusters contain all of the mass. How- 
ever, when t > 1, there is a second, nontrivial solution, 
< g < 1, and this solution is the physical one. Finite 
rings account for only a finite fraction, Mi = 1 — g, of the 
total mass. Therefore, giant rings must account for the 
remaining fraction of the total mass, g. At time t > 1, 
the total mass of the giant rings equals g{t)N . 

At time t = 1, the critical size distribution has a power- 
law tail (Figure , 



c fc (l) 



1 



k-">'\ 



(7) 



when k ^> 1. At the critical point, the size of largest ring 
scales as iV 2/3 with system size N [27l - [29| . 




FIG. 3: The critical size distribution, Cfe = Ck(t = 1) versus k. 
The simulation results are from 10 4 independent realizations 
of a system with N = fO 8 . 



IV. GIANT RINGS 




/ 

FIG. 4: Simulation results for G{1) = G(£, t) at three different 
times, t(g = 1/4) = 1.150729, t(g = 1/2) = 1.386294, and 
t(g = 3/4) = 1.848392. The time t(g) = ± In ^ follows from 
J6]). Also shown for reference is the theoretical prediction (|10l) . 
The simulations results are from 10 7 independent realizations 
of a system with size N = 10 6 . 



When t > 1, we expect that macroscopic rings, that is, 
rings that contain a finite fraction of the total mass in the 
system, account for the missing mass. For a giant ring 
with size k oc N, the total aggregation rate, k, and the 
total fragmentation rate, ojy ' are both proportional 
to N . Hence, aggregation and fragmentation occur with 
comparable rates. Also, since both rates are proportional 
to the system size, aggregation and fragmentation are 
very rapid. To find the size distribution of the giant rings, 
we must consider the aggregation-fragmentation process 
governing the giant rings. 

We characterize a giant ring using the normalized size, 
£, defined as £ — k/N . This quantity equals the fraction 
of the total mass contained in the ring. Let G(£, t) be the 
average number of rings with normalized size £ at time t. 
Conservation of mass dictates 



(8) 



g(t) = J dt£G{£,t), 

where g(t) is the nontrivial solution of ([6]). 
The quantity G(£, t) satisfies 



1 dG(£,t) _ 1 



N dt 



dss{£- s)G(s,t)G(£- s,t) 



- £(g~£)G(£,t) 

+ [ dssG(s,t)-^£ 2 G(£,t). (9) 

This rate equation, essentially the continuous analog of 
Eq. (j3|) , describes the aggregation- fragmentation process 
that governs the giant rings. To obtain © from ©, we 
first make the transformations Gk = Nrk and k = £N, 
and then, note that the aggregation loss rate is reduced 
by the factor (g — £) because self-interactions do not lead 
to aggregation. 



Since the total aggregation and fragmentation rates are 
both proportional to the system size, the left-hand side 
of (]9]) is negligible in the large- N limit. We therefore re- 
place the left-hand-side of Eq. ([9]) with zero to determine 
the time- dependent distribution G(£,t). The resulting 
non-linear integral equation has the remarkably simple 
solution (see Figured]) 



G(£,t) = 







£ < g(t), 
£>g(t). 



(10) 



Indeed, this solution obeys the mass conservation state- 
ment ((HJ). Surprisingly, the size distribution is universal, 
although the span of the distribution grows with time, 
< £ < g(t). Therefore, at time t > 1, there are giant 
rings of all sizes up to the maximal value g(t)N. 

The distribution of rings includes two distinct compo- 
nents: Nek gives the average number of finite rings, and 
G{£) gives the average number of giant rings. Of course, 
the former expression applies at all times, while the lat- 
ter holds only for t > 1. The giant rings grow at the 
expense of the finite rings and eventually, they take over 
the entire system as g — ¥ 1 when t —¥ oo. 

Finite rings and giant rings undergo separate, essen- 
tially decoupled, aggregation-fragmentation processes. 
Indeed, the rate equation (j4} for ct is closed, while the 
rate equation (|9|) for G{£) is, in practice, also a closed 
equation. There is a constant flux of mass, N x dg/dt, 
from finite rings to giant rings, and this flux couples the 
two types of rings. This coupling enters only through the 
fraction g(t) which appears explicitly in ([9]). 

The distribution (fT7j|) implies that there are multiple 
giant rings: the average number of giant rings, N g , scales 
logarithmically with system size (Figure [5J 



A^„ ~ In N. 



(11) 




N 

FIG. 5: The average number of giant rings, N g , and the vari- 
ance a 2 , versus system size N. The simulation results repre- 
sent an average over 10 5 independent realizations. We mea- 
sured N g and a 2 by counting the number of rings with size 
k > 41niV at time t = 2. 

This behavior follows from N g = ^dlG{t). The 
lower limit = k*/N can be deduced from the crite- 
rion NJ2k>k c k(t) = 1 that estimates the size of the 
largest finite ring. Using this criterion together with 
Eq. ([5]) we find fc* ~ (t — lni — In N and therefore 

Since the merger-breakup process is random, we ex- 
pect that the variance in the number of giant rings, a 2 , 
is proportional to the mean, a 2 ~ In AT. Numerical sim- 
ulations confirm this behavior (Figure [5]). Hence, the 
standard deviation 

a-Vh^N (12) 

quantifies fluctuations in the number of giant rings. 

Figure [B] shows the normalized sizes of the three largest 
rings as a function of time using a simulated system. 
These sizes exhibit huge fluctuations as giant rings con- 
stantly merge and break on a very fast time scale. In- 
terestingly, while the size of an individual giant ring is a 
stochastic quantity, the total size of all giant rings, g(t), 
is a deterministic quantity. 

The number of finite rings is proportional to N, while 
the number of giant rings scales only logarithmically with 
N. Equation (j^J shows that monomers dominate in the 
long-time limit. By comparing the average number of 
monomers, Nci = Ne'*, with the the average number of 
giant rings, given by we conclude that giant rings 

overtake finite rings when t 3> tf with 

tf~1nN. (13) 

In writing this expression, we ignored secondary logarith- 
mic corrections. 

For times t 3> tf, the ring-size distribution reaches a 
steady state. Setting g = 1 in (fTOj) shows that Nk, the 
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FIG. 6: The largest three rings. Shown is the time evolution of 
tn{t), the size of the nth largest ring at time t, for n = 1, 2, 3. 
The results are from a single run of a system with N = 10 6 . 
Also shown is the cumulative mass git). 

average number of rings with (unnormalized) size k, has 
the following form 

N k = i, (14) 

for all 1 < k < N. Thus, at the steady state, there are 
rings of all lengths, from finite rings to macroscopic rings. 
The distribution (|T4")) also follows from the detailed- 
balance condition KijCiCj = FijCi+j with the solution 
c k = (Nk)- 1 [Ull. 

V. SHUFFLING ALGORITHM 

Throughout this paper, we presented results of Monte 
Carlo simulations that support the theoretical predic- 
tions. We implemented an elegant algorithm that takes 
advantage of an isomorphism between the polymerization 
process and a card shuffling process. In the card shuffling 
algorithm [3lTj34l |. we start with an ordered deck of N 
cards, numbered 1 through N. Then, at each elementary 
step, we pick two cards at random and exchange their 
positions. For example, the first two step in shuffling a 
deck of 6 cards may look like 

123456 -> 153426 -> 154326 

Time is updated by At = 2/N after each step, 
t —> t + 2/N, and thus, each card participates in one shuf- 
fling event per unit time, on average. 

We now use cycles to represent permutations. For 
example, six cards ordered 134265 are represented by 
(1) (234) (56) because there are three cycles: the card 1 
forms a cycle of length one, the cards 234 form a cycle of 
length three, and the cards 56 form a cycle of length two. 
Initially, there are cycles of length 1. Then, exchange 
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FIG. 7: The average number Nk of rings of size k at the 
steady-state. The simulation results are from 10 3 independent 
realizations of a system with N = 10 8 at time t = 20. 

of two cards in distinct cycles leads to merger, while ex- 
change of two cards in the same cycle leads to breakup. 
For example, the following steps generate the merger and 
breakup events in Figure [2j 

(123) (456) (156423), and (156423) (123)(456). 

Furthermore, the merger rate and the breakup rate are 
given by ((XJ) and ©. Hence, the dynamics of cycles in the 
shuffling process arc identical to the dynamics of rings in 
the polymerization process. 

The above algorithm is straightforward and efficient. 
The shuffling steps take O(N) operations per unit time, 
and moreover, tracing the cycle structure requires only 
O(N) operations. This linear algorithm enabled us to 
simulate large systems with as many as N = 10 8 . Fig- 
ure [7] demonstrates the excellent agreement between the 
simulation results and the theoretical prediction ([T4|. 

The distribution Nk given in (fT4|) equals the average 
number of cycles of length k for a random permutation of 
N elements [35|,|36j]. As expected, repeated shuffling ran- 
domizes the card order and according to (|13l) . the num- 
ber of exchanges required to generate a perfectly random 
shuffle scales as N In N. 

A natural generalization is to n-card shuffling where n 
randomly chosen cards are reordered according to a pre- 
scribed rule. For example, if n — 3, we may follow the 
cyclic rule 123 — > 231. The equivalent polymerization 
process now involves merger of n polymers. Straightfor- 
ward generalization of the F lory model shows that the 
total gel mass, g(t), satisfies [37| 

1-5 = e ^Ti— *, (15) 



with < g < 1 in the giant- ring phase t > (n — 2)1. We 
anticipate that the distribution of giant rings is given 
by (fl"0|) . except that the total mass is now specified by 
(|15p . Our numerical simulations of the three-card process 
confirm this behavior. 

VI. DISCUSSION 

In summary, we studied a ring polymerization process 
in which a bond between two monomers results either 
in aggregation of two rings into one or in fragmentation 
of one ring into two. This process exhibits a percola- 
tion transition with a finite-ring phase in which all rings 
are microscopic and a giant-ring phase in which multiple 
macroscopic rings coexist. While the cumulative mass of 
the giant rings is deterministic, the sizes of individual gi- 
ant rings are stochastic. Moreover, the giant rings exhibit 
huge fluctuations due to the extremely rapid merger and 
breakup processes. Finally, the size distribution of giant 
rings is stationary, although the span of this distribution 
grows with time. 

The aggregation-fragmentation process that governs 
the rings is perfectly reversible. On the one hand, the 
distribution of ring size reaches a stationary state where 
detailed balance is formally satisfied. On the other hand, 
this final distribution is not thermodynamic because the 
number of rings varies logarithmically, rather than lin- 
early, with system size. Phase transitions with non- 
thermodynamic states were previously observed only in 
irreversible aggregation- fragmentation processes [38l - l40| . 

The ring formation process can be generalized in many 
ways. We focused on the mean-field version, and it will be 
interesting to study two-dimensional rings where spatial 
correlations play an important role. Another direction 
for further study is percolation of polymers with other 
types of fixed structures, for example, polymers where 
all monomers have exactly three bonds [3j. 

Finally, we notice that the unusual behaviors in the 
giant-ring phase are the consequence of the basic topo- 
logical constraint, namely that the polymers maintain a 
ring structure. This suggests to investigate the influence 
of other constraints such as planarity [4ll . [42| . Another 
interesting ques tion is what happens if the polymers are 
membranes [43[ such as spheres (say due to the surface 
tension) that can merge and divide. 

We thank Kipton Barros for useful discussions and 
Talia Bcn-Naim for experimenting with magnetic rings. 
This research is supported by DOE grant DE-AC52- 
06NA25396 and NSF grant CCF-0829541. 
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